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, Abstract. We consider the steady state of a one dimensional diffusive system, such as the sym- 

metric simple exclusion process (SSEP) on a ring, driven by a battery at the origin or by a smoothly 
varying field along the ring. The battery appears as the limiting case of a smoothly varying field, 
^ . when the field becomes a delta function at the origin. We find that in the scaling limit, the long 

^ ' range pair correlation functions of the system driven by a battery turn out to be very different from 

, the ones known in the steady state of the SSEP maintained out of equilibrium by contact with two 

I ■ reservoirs, even when the steady state density profiles are identical in both models. 
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' There are notoriously few fully analyzable models of non trivial (interacting) current carrying 

systems (Lj [Sc] in non-equilibrium steady states (NESS). The few exceptions are almost all one 
dimensional lattice gases evolving according to stochatic jump processes and interacting via exclu- 
^ ■ sions, in contact with particle reservoirs at different densities. Among the simplest of such models 
O . is one of particles on a lattice of N sites (i = 1, . . . , N) in which the bulk dynamics are given by the 
symmetric simple exclusion process (SSEP) and sites 1 and N are in contact with particle reser- 
voirs. In the case of reservoirs at equal densities the stationary state satisfies detailed balance. For 
reservoirs at unequal densities, say pa on the left and pb on the right {pa > Pb), the matrix method 
QQ ' gives a full microscopic description of the NESS including the explicit form of the correlation func- 
. tions Pkih, ■ ■ ■ ,ik', Pa, Pb, |DLS2| and allows one to calculate the large deviation function (LDF) 
0^ ; in the hydrodynamical scahng hmit, N ^ oo, i/N ^ x e [0,1] |DLSH IBDO.TLII |BDGJL2]. One 
can also obtain, in this scaling limit (with diffusively scaled time), the time evolution of the system 
! in a non stationary state via a diffusion equation for typical density profiles as well as the LDF for 
^ ; time evolving densities and currents |BDO,TLll IBDG.TLaj iPTllBm] . 

In attempts to further extend our understanding of such non-equilibrium systems, we investigate 
here the stationary state of a system of M particles on a ring of sites with exclusion. The 
^ . jump rates to the right and left (counter clockwise, clockwise) for a particle at site i are pi and qi, 
^1 i = 1 . . . N. We shall be interested, as usual, in the case N ^ 1, M/N = p. We consider two types 
of situations: 

(1) The "battery" case: the system evolves according to the SSEP except that there is a 
"battery" at bond (A^, 1), i.e. = 1 for i / A^, = 1 for z / 1 and pN = p,qi = q with p 
and q independent of N. 

(2) The smooth asymmetric case (WASEP) Pi = 1 + jjE{i/N) and qi = I - jjE{i/N) with E 
a smooth field of period 1. 

(Combination of both cases can also be treated). 
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In the first case, the microscopic "battery" at (A'", 1) induces a NESS with a particle current of 
magnitude proportional to 1/N and a concomitant linear density profile. This looks very similar 
to what happens in the open system and would suggest that the measures for the reservoirs and 
battery driven systems would be similar when l:^ 1, in the same spirit as the difference between 
"canonical" and "grand-canonical" nature of the two systems disappears in equilibrium, when p = q 
and M/N Pa = Pb- This is indeed the case to the leading order in N. Local equilibrium holds in 
both models and implies that both steady state statistics are given locally by a product Bernoulli 
measure. Here, we will focus on the corrections to this local Bernoulli measure and on the long range 
correlations. We will see that the two NESS measures differ much more than the corresponding 
equilibrium ones. The complicated interplay between the non-equilibrium long-range correlations, 
the canonical constraint and the driving mechanism leads to two-point correlation functions with 
different structure for the reservoir and battery driven systems. In particular, the long range 
correlations are much more singular for the battery model than those in the open system. A 
variant of this model with a macroscopic battery was introduced in |IIF| and a related dynamics 
with colliding hard spheres was also considered in |BPSj . 

In the second case, the macroscopic stationary density profile p{x) then satisfies the general 
equation for driven diffusive systems 

d,p{x) - 2E{x)p{x) (1 - p{x)) =-J, (1.1) 

where J is the stationary current which is independent of x. Equation (jl.ip corresponds to the 
stationary solution of the diffusion equation with a drift 

dtp = dip-dAm^)p{^-p)) ■ (1-2) 

The first case can be thought of as a very singular limit of the second one when E[x) becomes a 
5-function localized at the origin. 

We will also consider the Zero Range Process (ZRP) driven by a "battery". This microscopic 
dynamics is simpler than the SSEP driven by a battery, as the NESS of the ZRP on the ring is 
a product measure constrained to having a fixed number of particles. There are no long range 
correlations in the non-equilibrium measure of the ZRP with reservoirs |MF] and the long range 
correlations of the ZRP driven by a battery are similar to those found for an equilibrium system 
with the constraint that the total particle number is fixed. 

The outline of the rest of the paper is as follows. In section [21 we define the microscopic models. 
By using a macroscopic approach, the two-point correlation functions are computed for the battery 
model in section [3] and for a slowly varying field in section [H These results are then compared, 
in section O to the two-point correlation functions of an open system in contact with reservoirs. 
Finally the invariant measure of the ZRP driven by a battery is computed in section [6l Some 
technical details concerning the solution of the macroscopic equations are given in the appendices. 

2. The models 

2.1. The battery model. We consider the SSEP on the ring {1, A^} containing M particles with 
jump rates 1 to the left and to the right except at the bond (A^, 1) where the jump rate from A^ 
to 1 is p and from 1 to A^ is q. These modified asymmetric rates act as a battery which forces a 
current of particles through the system. We assume that p > and q> 

In the absence of the battery [p = q = 1), the stationary state is one in which all the (^a/^ 
configurations have equal weight. For A^ — )• cxd, M = pN, this corresponds locally to a product 
Bernoulli measure with density p. Denoting by the Bernoulli measure at site i with density p, 
the product measure (^jg^ invariant wrt the SSEP dynamics for any constant density p. 
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For the SSEP on Z with a battery at the bond (0, 1), the invariant measure remains a product 
but the density is discontinuous across the battery. The invariant measures on Z are 



^^"'^" = (g)^r (8)'^r, (2.1) 

j ^ i > 1 

where the densities p~ , p'^ (at sites and 1) satisfy the equation 

pp-{l-p+)=qp+{l-p-). (2.2) 

These measures satisfy detailed balance and there is no current flowing in the system. 

On the ring when p ^ q, the steady state measure is not known. Taking averages with respect 
to the time evolving measure ij,n,t, we obtain 

ij^l,N, drPN,T{'ni) = AiAf,r(??i+l) + /iAr,r (??i-l) " 2pN,r{'ni) (2-3) 

drPN,T{'ni) = PN,T{m) - PN^rim) + ppN,T{'nNii - m)) - QPN^rivii'^ - vn)) 

dTPN,T{'nN) = Pn,t{vn-i) - fJ'N,T{r]N) - PPN,T{riN{^ -Vi)) + 9/"A',t(??i(1 - Vn)) 
It follows immediately that in the stationary state, the average density profile is linear 

im) = im) + ^;f^{{m) - im)) , (2.4) 

where (•) stands for the stationary measure and the mean density p = ^ satisfies 

P=1{{vn) + {vi))- (2.5) 



There is also a stationary average current q of order 1/iV 

1 



^{{m)-{VN)), (2.6) 



and one has the identity 

p{vn{1-Vi)) -QiviC^-m)) =Q- (2.7) 
For finite A^, the relations (j2.4p and (|2.7p do not allow to determine the average profile (r/j) as (j2.7p 
involves the correlation between sites N and 1. If one assumes however that local equilibrium holds 
for large then the system will behave at the battery as if it was described by the measure ()2.ip . 
Thus one expects that for large 

{riNVi) ~ iVi) (Vn) = P'^P~ , (2.8) 

with p^ satisfying 

pp-{l-p+)=qp+{l-p-), P^±^ = -p. (2.9) 

where we used (|2.5p and (|2.7p . This implies that in terms of the macroscopic variable x = i/N, 
that the steady state density p{x) is linear with a discontinuity at the battery 

p{x) = p+ + {p- - p+)x, with p{0) = p+^ p{l)=p-. (2.10) 

The local equilibrium (12. Sp can be justified [BL] and one can prove that after space/time rescaling, 
the microscopic equations lead to a macroscopic description of the local density by the heat equation. 
The battery yields non-linear boundary conditions 

' dtp{t,x) = dlp{t,x) 

p p{t,l){l - p{t,Q)) =gp(t,0)(l-/>(t,l)) (2.11) 
[a^.p(t,0) =a^/j(t,l) 



Vt > o,x e]o,i[, 
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where p{t, x) stands for the local density at the macroscopic time t and position x. 

2.2. The case of a slowly varying field. In contrast to the battery where the driving force is 
localized on a single bond, we will now consider the case of a smoothly varying macroscopic field 
E{x) along the ring. The microscopic model is now the weakly asymmetric simple exclusion process 
of M particles on a ring of sites. A particle on site i jumps to its neighboring site on its right at 
rate pi and to its neighboring site on its left at rate qi provided that the target site is empty. By 
weak asymmetry, we mean that pi and qi have the following scaling dependence on the system size 
N 

When the integral of E{x) over the circle is non zero, the system will reach a non-equilibrium 
steady state for every fixed A'^. 

The time evolution of the microscopic mean density is given by 
<9T/iAr,T(??j) = Pi-itiN.T{vi~i{'^-Vi))-QifJ'N,T{r]i{l-r]i-i))-piPN,T{m{'^-'ni+i))+Q^^^ 

(2.12) 

where /iAr,r(') is the expectation value wrt the time evolving measure. 

Assuming that for large A^, the density profile takes the following scaling form 



and the correlation functions scale as 



1 I ^ ^ J \ 



2 



one can show that for x = i/N, t = t/N 

Pi-iiJLN,T{'ni-i{'^ - m)) - qifj'N,T{'nii^ - 'm-i)) 

~ ^ [-d,p{t, x) + 2E{x)p{t, x){l - p{t, x))] + -^Wit, x) , 

where W depends on the two-point correlation function C defined in ()2.14p . One then gets from 
(j2.12p the macroscopic evolution equation 

dtp{t,x) = d,[d,p{t,x) -2E{x)pit,x)il - pit,x))] . (2.15) 

This is the viscous Burgers equation on a ring x £ [0, 1]. A mathematical derivation of (j2.15p can 



be found in \KL\ Sp2 



The macroscopic evolution of the battery model (j2.1ip can then be recovered by taking a large 
localized field E{x) at the origin with a given integral equal to K. According to (|2.15p this leads, 
in the limit, to a jump of density across the origin 

1-p Jp- p(l-p) 

This expression is equivalent to (j2.2p with 2K = log{p/q). The effect of the battery should be 
understood as a force which maintains a fixed difference of chemical potentials across the origin. 
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2.3. The Zero Range Process and the moving battery model. In section EIJ we will study 
the invariant measure of a ZRP driven by a battery. The microscopic dynamics is defined as follows. 
At site i, the occupation number r/j can take any integer value and a particle at site i performs a 
jump to the left or to the right with rate g{r]i), where the function g is increasing and g{0) = 0. 
The driving force is modeled by modifying the jump rates between sites 1 and N: a particle jumps 
from 1 to with rate q g{i]i) and from A to 1 with rate p girjN)- The hydrodynamic limit for this 
model on Z was derived in |LOVj . 

The ZRP driven by a battery is related to the SSEP on a ring driven by a "moving battery", 
i.e. a special tagged particle with asymmetric jump rates p and q. One can think of this particle 
as being driven by an external field E with p/q = exp(2E) |FGLt [LOV| . A current will be induced 
by the moving battery. The mapping between the two models was exploited by jLQVj to study the 
motion of this tagged particle in Z. More precisely, the SSEP with a moving battery and with M 
particles on the ring of length N can be mapped onto a ZRP on a ring of length M with N — M 
particles and with the specific rates g{n) = In ^ i- Let 1 be the index of the asymmetric particle in 
the SSEP and 2, . . . , M the indices of the other particles, then the ZRP variable rji stands for the 
number of empty spaces ahead of particle i. 

3. The macroscopic approach for the battery model 

A generic property of non-equilibrium systems, maintained in a steady state by contact with 
reservoirs at unequal chemical potentials, is the presence of long range correlations |Spl[ISC[IDKSl 
IDELOl [BDGJL4]. The steady state correlations can be predicted from a macroscopic approach 
based on " fluctuating hydrodynamics" |Spl[ IDELOj IPS] . Alternatively, one can derive the density 
large deviation functional for the steady state and then recover the correlations by expanding 
the functional near the steady state [BDGJLll IBDGJL4t IDT| IBDLWj . We will apply the latter 
approach to compute the correlations in the battery model introduced in section 12.11 

3.1. The hydrodynamic large deviations. After rescaling the space hy N (i = xN) and the 

time by A^ (r = NH), the microscopic system can be described at the hydrodynamic scale by two 
macroscopic functions the local density p{t, x) and the local current of particles q{t, x) which obey 
the conservation law 

dtp{t,x) = -d^q{t,x) . (3.1) 

The probability of observing a joint deviation of the current and the density in our microscopic 
system depends exponentially on the system size. More precisely, the probability of observing an 
atypical macroscopic trajectory {p{t, x),q{t, x))o sg * ^ t during the macroscopic time interval [0, T] 
scales like 

^[o,TN'^](^{p{t,x),qit,x))^ ~ exp (^-AJ"[o,r](p,g)) , 
where the large deviation functional is given by 

rT rl (^{t,x)+d^p{t,x)f 



J-[0,., (p, ,) = dt dx^^^^ . (3.2) 

Here a{u) = 2u{l — u) and the density has to satisfy the conditions 

pp{t,l){l-p{t,0))=qp{t,0){l-p{t,l)), f dxp{t,x)=p. (3.3) 







If {p{t,x),q(t,x)) does not satisfy (j3.ip or if the conditions (j3.3|) are not satisfied (for a set of 
times with non-zero Lebesgue measure) then the functional is infinite. The expression (13. 2p is a 
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generalization of the large deviation functionals derived for open systems (see |BDGJL2| IDl^ IBD2] 
for reviews). The condition ()3.3p at the battery follows from (I2.2p . as local equilibrium is still 
satisfied in the hydrodynamic large deviation regime. 

To compute the probability of observing an atypical density trajectory x))o ^ t ^ t 

IP[o,TAf2] (^Pit, x)^ - exp (-iV-F[o,T] (p)) , with J'p^T] (p) = inf -^[o,T] (p, o) , (3-4) 

one has to optimize J^[o,t]{P:Q) over all the currents which are compatible with the density (|3.ip 
and which can be written as 

rx 

Qit,x)=j{t)- dtp{t,u)du, (3.5) 
Jo 

where j{t) is the current at 0. Optimizing ()3.2p over j{t) implies 

kp{t, u) du 
a{p{t,x)) 



We introduce an auxiliary function H such that 

q{t, x) = -d^p{t, x) + a{p{t, x))d^H{t, x) , (3.7) 

and get from 



0= r = [\xd,H{t,x) = H{t,l)-H{t,0). 

Jo (^{p{t,x)) Jo 

This leads to a continuity condition on H at the battery 

H{t,0) = H{t,l)- (3.8) 

For a trajectory (p(t, x))o ^ t ^ t, there exists H satisfying (j3.7p and such that 

dtp{t, x) = dlp{t, x) - (a{p{t, x))d^H{t, x)) , (3.9) 
thus the functional (13.41) reads 



J'[o,T]{p) = ^ I dt [ dxa{p{t,x)){d^H{t,x)f . (3.10) 







2 

Note that finding the optimal trajectories {H,p) is equivalent to finding the optimal {q,p)- 



3.2. Steady state large deviations. Following [BDG JLT] . the density large deviations for the 
steady state can be computed by using the hydrodynamic large deviations. For any smooth function 
A(x) in [0,1], one defines 

which is given according to |BDGJLT] by 

g{X) = lim sup I / A(x) p(0, x) - Ti^T,o] (p) | , (3.12) 

where J~[-t,o] is defined by (j3.10p . The supremum ranges over all the density trajectories p{t,x) in 
the macroscopic time interval [—T, 0] which are equal to the steady state p{x) at time —T. 
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Given the function A and a time — T, we look for the optimal density trajectory p in the variational 
problem ()3.12p . To determine this optimal trajectory, we consider a variation p ^ p -\- (p and 
H H + h. From (|3.9|) . <p and h satisfy 

dtv = dl^ - {a'{p)H' ^ + a{p)h') . (3.13) 



To simplify notations, we have omitted the (t, x) dependence and used the shorthand H' = 
The conditions (I3.3|) become 

^^'''^ f\x^{t,x)=0. (3.14) 



p{t,i){i-p{t,i)) p(t,o)(i-/>(t,o))' 7o 

As the density is equal to p at time — T then ip(—T, x) = 0. 

For an optimal trajectory p, the first order perturbation in ip in the variational problem ()3.12|) 
should be equal to so that 



^ dxX{x)ip{0,x) - J dt dx (^^a'{p){H'Y ip + a{p)H' h'^ =0. 



(3.15) 



Integrating by parts and using (j3.13p . one gets at time t 
dx a{p)H'h' 



1 

dx {d,{a{p)h'))H + {Hit,l)a{p{t,l))h'it,l) - Hit,0)aipit,0))h'it,0)) 

= dx (^dt^ - dl^p + a, (a'(p)Fv)) a 

J 

+ (i/(t, l)a(p(t, l))/i'(t, 1) - Kit, 0)a(/5(t, 0))/i'(t, 0)) . 
Another integration by parts leads to 

dt dxa{p)H'h' = - dt dx [ipdtH - <p' H' + ipa' {p){H'f] (3.16) 

J-T Jo J-T Jo 



+ [ dx ip{0,x)H{0,x) + [ dtH{t,l)j{t,l) - H 

Jo J-T 



(t,0)i(t,0) 



where j(t, x) = —(p'{t, x) + cr'{p{t, x))H'{t, x)ip{t, x) + a{p{t, x))h'{t, x) is the current variation q — )■ 
q + j. The current j{t,x) is continuous at the battery (this can be checked by combining the 
condition dxtp{t,x) = and (|3.13p which says that dt(p = —j')- As H is also continuous at the 
battery (j3.8p . the last boundary term in (j3.16p vanishes. One has 

t-O r-l rO rl 

dt dxa{p)H'h'= dt dx [p {-dtH - dlH - a' {p){H'f)] 

J-T Jo J-T Jo 

+ I dx ip{0,x)H{0,x) + [ dt [H'{t,l)p{t,l) - H'{t,0)ip{t,0)] . 

Jo J-T 
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Thus for any perturbation (p, the condition ()3.15p can be rewritten 

fO /-l 



dt / dx 

T 



^[dtH + dlH + ]^a'{p){H'f 



+ / dx [\{x) - H{0, x)) v?(0, x) 
Jo 

- dt [a{p{t, l))H\t, 1) - a{p{t, 0))H'{t, 0)] = , (3.17) 

where we used the identity (j3.14p . 

Combining (|3.17p and (j3.9p . the evolution equations for the optimal trajectory in (|3.12p are 

\m = -d'^H-la'ip){H'f 

with boundary conditions 

Va;G[0, 1], p{-T,x) = p{x), H{0,x) = \{x). (3.19) 
We also have that at the battery, for any time t < 0, 

pp{t,l){l-p{t,0))=qp{t,0){l-p{t,l)), p'{t,0)=p'{t,l), (3.20) 

H{t, 0) = H{t, 1), p{t, 1)(1 - p{t, l))H'{t, 1) = p{t, 0)(1 - p{t, 0))H'{t, 0) . (3.21) 

The boundary conditions fohow from ()3.3p . (j3.8p and (j3.17p . The continuity of p' in (j3.20p is a 
consequence of the continuity of the current (j3.7p at the battery and of (j3.2ip . We stress the fact 
that H{0,x) = A(x) in (|3.19p may not satisfy the boundary conditions ()3.2ip . 

The bulk evolution equations (j3.18p have already been derived in previous works [BDGJLTj 
IBDLWl iDGj . but the boundary conditions (j3.20p . (j3.2ip are specific to the battery model. 

3.3. Linearized equations. The two-point correlations in the steady state can be obtained by 
taking the second derivative of ^(A) at A = (see (|3.1ip ) |BDGJLl| IDlj . This can be understood 
as follows. For a given A^, one has at the second order in A 

(3.22) 

' iv^Kl - {v^)?j + ^ E A (^) A (^) ivuV,) + o(A2) , 

where {r]i]r]j) = {rjiVj) ~ iVi) iVj) is the connected two-point correlation function. The above 
expansion amounts to perturbing the system around the steady state p{x) (|2.10p . 
For N large, one expects from (j3.1ip and (j3.22p that 

g{X) = C dx -p{x)\{x) + -p[x)(l - p(x))^ + /' dx C dy C{x, y)X{x)X{y) , (3.23) 

Jo ^ JO Jx 

where C{x,y) is the macroscopic non-equilibrium two-point correlation function 

at the leading order in A^. C{x,y) will be computed in (j3.37p . We stress the fact that the equality 
()3.23p rests on the assumption that one can interchange the limits N ^ oo and A ^ 0. 

Our goal now is to compute ^(A) to the second order in A from the variational problem (j3.12p and 
to deduce from it C{x, y) (j3.23p . To any A, one can associate an optimal density px which minimizes 
the variational problem (I3.12p . By construction the density at time satisfies px{0,x) = so 
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that expanding px{0,x) to the first order in A will lead to the second derivative of Q{X)- One can 
see from (I3.23P that 

px{0,x) = p{x) + p{x){l-p{x))X{x)+ [ dyC{x,y)X{y). (3.25) 



Jo 

To determine p\ for small A, we linearize the optimal evolution ()3.18p around the steady state 
profile p (|2.1U|) . For the sake of notation, we drop the subscript A in px and decompose the trajectory 
as p\ = p + (f and H = h, where ip and h are small. At the first order, ()3.18p becomes 

vt€[-r,o], ldt^ = d'.^-d.{a{p)h') ^3^^^ 

[dth = -d'^h 

and the boundary conditions (13.190 . (I3.20p . (I3.2ip lead to 

VxG[0,1], ip{-T,x) = 0, h{0,x) = X{x), (3.27) 

and 



a{p ) cj(p+)' 



ip'it,0)=ip'it,l), (3.28) 



h{t,l) = h{t,0), a{p-) h'{t,l) = a{p+) h'{t,0) . (3.29) 

Note that ()3.28p has been obtained by linearizing ()3.20p and that h is defined up to a constant 
which does not influence the evolution of the density. 

The coupled equations (|3.26p can be solved by integrating the heat equation. We first start 
with the second equation. The Green's function associated to the Laplacian with boundary 
conditions (j3.29p satisfies for any y in (0, 1) 

dtGUx,y) = d',GUx,y), with G?(l, y) = G?(0, y), aa,G?(l, y) = a,G?(0, y), (3.30) 

where we set 



aip+) 



(3.31) 



The function h follows the backward heat equation on [— T, 0] with final data /i(0, y) = X{y) at time 
0. Therefore, we can write 



VtG[-T,0], h{t,x)= [ Glt{x,y)X{y)dy. 

Jo 



We turn now to the evolution of in (|3.26p which follows the heat equation with a source term 
depending on h. The Green's function G^ associated to the Laplacian with boundary conditions 
satisfies for y in (0, 1) 

dtGjix,y)=d^,Gjix,y), with Gjihy) = aGjiO^y), d,Glil,y) = d,GjiO,y) . (3.32) 
Given Lp{—T, y) = at time — T 



Vt € [-r,0], ^{t,x) = - j ^ds dy GU{x,y)dy(a{p{y))dyh{s,y)) 



10 



T. BODINEAU(i), B. DERRIDAt^), J. L. LEBOWITZ^^) 



One can check directly that the Green's function solution of (I3.32p is given by 



2 4 °° / 
Gl{x,y) = -a{x) + -^exp {- {2kTT ft) (a{x) cos{2tt kx) cos{2tt ky) + (3.33) 

k=l 

a{l — y) sin(27r/c2;) sin(27r/c?/) + 4/i:7r(l — a)t sin(27r/cx) cos(27rA;?/)^ , 

with a{x) = ax + \ — X and a = O^i^ can also check that the Green's function which 

solves (j3.30p is given by 

G]{x,y) = G][y,x). (3.34) 
In appendix I, the construction of the Green's functions is explained. 



We are going to relate (^(0, x) and A(x). The density fluctuation at the final time reads 
(/?(0,x) = - I ^ds j dy j dwG\^{x,y)dy{a{p{y))dyG'^_s{y,w)^\{w), 



which is equivalent to 



99(0, x) 



dw 



ds / dyGl{x,y)dy[a{p{y))dyGl{y,w] 
Jo ^ 



X{w). 



In appendix II, it is shown that (j3.35|) can be rewritten as 



1 



ip{0,x) = -a{p{x))X{x) + / dyG{x,y)X{y) 



where the long range correlations are given by 

2 



Gix,y) 



a(p+)+a(p-) , 
, ^' ' ^'^ ' + — ]{ax + l-x)iay + l-y) 
(a + 1)2 V 2 3 7^ ^ ^' 

-2J^ ds dzGlix,z)Gliy,z)-Gl{x,0)Gliy,0). 
Jo Jo 



Since px{0,x) = p{x) + 93(0, x), this determines the correlation function in (j3.25p . 



(3.35) 



(3.36) 



(3.37) 



4. The macroscopic approach for the slowly varying field 

In section [2?2l a microscopic dynamics driven by a weak field E{x) was introduced. We consider 
now a general diffusive system at a macroscopic level and write down the two-point correlation 
function which generalizes p.37p to an arbitrary weak field E[x). In the limit where E{x) becomes 
a (5-function, we will see that one can recover the correlations of the battery model ()3.37p . 

The main advantage of the calculations of this section is that the effect of the battery is 
smoothened over the whole system through the function E{x), so that there are no longer bound- 
ary conditions such as (j3.20p or (j3.2ip at the battery and the functions H{t, x) and p{t, x) become 
simply periodic functions of the space variable x. One difficulty is that for a generic E(x) the 
steady state profile ^(x) and the Green's functions G^ and G^ are not known explicitly and we 
have to make a few assumptions on the convergence of hydrodynamic equation (14. 5p to the steady 
state p{x) or on the long time behavior of the Green's functions (I4.18p . 
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4.1. The variational problem. Our starting point is that the large deviation functional (|3.2|) of a 
macrosopic trajectory {p{t, x), q{t, x))q ^ t ^ t which, for a general diffusive system is characterized 
by two functions D(p) and a{p), takes the following form [BDGJL3|, IBD2j in presence of a driving 
field E 

fT ri (q{t,x) + D{p{t,x))d^p{t,x)-E{x)a{p{t,x))y 
•^[o,T, (P, <!) = I dt dx ^^^^^ , 

The SSEP with a weak field introduced in section [2^2] corresponds to the functions D{p) = 1 and 
a{p) = 2p{l — p). Following exactly the same steps as in subsection 3.1, one gets, by optimizing 
over the current defined in (|3.5p . that ()3.7p becomes 

q{t, x) = -D{p{t, x)) d^p{t, x) + a{p{t, x))[E{x) + d^H{t, x)] , (4.1) 

where H(t,x) is periodic in space {H(t,x) = H(t,x + 1)). This implies that the time dependent 
density profile p is related to H by 

dtpit, x) = {D{p{t, x)) d,p{t, x)) - x)) [E{x) + d,H{t, x)]) , (4.2) 

(instead of (|3.9p ) and the density large deviation functional is given by 

•^[0,T](P) = \ dt j\x a{p{t,x)){d.H{t,x))\ (4.3) 

as in (f3lB . 

To evaluate Q defined in (j3.1ip . (j3.12p we proceed as in subsection 13.21 and determine the optimal 
dynamical fluctuation starting from the steady state p{x) given as the solution of 

- D{p{x))-p'{x) + E{x)cj{-p{x)) = J , (4.4) 

where JT" is the steady state current. In Appendix III, we check that for regular coefficients 
a{p),D{p) > and E{x) considered here (j4.4p has a unique solution. We shall further assume 
that the hydrodynamic evolution 



dtp = [D{p)d.,p - E{x)<j{p) ] , (4.5) 
converges to the steady state profile p{x) 



As in subsection 13.21 one considers a small variation p^p + (p,H^H + hoi p{t, x) and of 
H(t,x) and one gets from (j4.2p that (|3.13p becomes 

dt^ = dl^Dip)^) -d,[E{x)<y'{p)^ + a'{p)H'^ + a{p)h!) . (4.6) 

Then following exactly the same steps as in subsection 13.21 one ends up with ()3.18p replaced by 

dtp = (D{p)d^p - E{x)a{p) - a{p)d^H] (4.7) 

1 

— ( 
2 



dtH = -D{p)dlH - E{x)cj'{p)d^H - \a'{p){d^Hf (4.8) 



with the boundary conditions at times —T and as in (j3.19p 

Va;G[0,l], p{-T,x) = p[x), H{0,x)=X{x). (4.9) 

The spatial boundary conditions (|3.2Up . (|3.2ip are now replaced by the requirement that H{t,x) 
and p{t,x) are smooth periodic functions of the space variable x (except possibly at time t = 0). 
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4.2. Small variations of the density. In order to determine the correlation function CE{x,y) 
as in ()3.25p we need to solve the above equations ()4.7l - lir9l) to first order in A(x), that is to first 
order in if = p — p and h = H. Linearizing ()4.7p . ()4.8p . one gets 



dt^ = d,\d,{D{p)ip) - E{x)a'{p)ip - a{p)d,h) (4.10) 
dth = -D{p)dlh - E{x)(t{p)' d^h (4.11) 

As in section [31 these equations can be solved in two steps. As the equation for h does not involve 
if, it is determined in terms of the Green's function as 



VtG[-r,0], h{t,x) 



G\{x,y)\{y)dy. 



where is solution of 



dtG^^{x,y) = D{-p{x)) dlGl{x,y) + E{x)a' {-p{x)) d^Gl{x,y) 
with Gq(x, y) = 6x=y To solve for ip, it is convenient to write f as 



(p{t,x) 



a{p{x)) 



2D{-p{x)) 

The evolution equation of can be determined from (14.8 



h{t, x) + ^(t, x) . 



(4.12) 
(4.13) 
(4.14) 



a'p' 



Eaa' 
2D 



where the functions D,a,a' are evaluated at density p{x) and E at position x. Using the fact that 
the steady state profile p{x) satisfies (I4.4p . this becomes 

dti^ = dl [Dij] - dx [Ea'ij] - Jhdx 



a 



2D 



We introduce the Green's function G\ as the solution of 

dtG\{x,y) = dl{D{-p{x)) G\{x,y)) - dx{E{x)a' {-p{x))G\{x,y) 
with Gq{x, y) = 6x=y Then one can solve (|4.15p 



(4.15) 



(4.16) 



(pix,0) 



2D{-p{x)y 



X{x)-J r dt I dz I dy G\ {x, z)d, ( ) Gl (z, y)X{y) (4.17) 



2D{-p{z)). 

dz I dy G\{x,z)^^^^G\{z,y)\{y) . 



We assume for all x the convergence when T oo 

G]^{x,z) = Gl{z,x) ^ B{x). 

For the battery model B{x) has the explicit form ^pj^ (see p.SSp ). Taking the limit T 
()4.17p gives the correlation functions 



(4.18) 
> oo in 



CE{x,y) = -B{x)B{y) / dz 



2D{-p{z)) 



J 



dt I dz G\{x, z)dz 



2D{-p{z)) 



GKz^y). 



As for the battery model, the Green's functions satisfy at any time t the symmetry property (|3.34p 
(see the proof below) 

Gl{x,y) = G]{y,x). (4.19) 



A DIFFUSIVE SYSTEM DRIVEN BY A BATTERY OR BY A SMOOTHLY VARYING FIELD 



13 



Thus we finally get the macroscopic expression for the two point correlation functions with a weak 
field 

CE{x,y) = -B{x)B{y) j dz^^^^ - J dt j dz G\{x,z)G]{y,z)d-\^^^) . (4.20) 



We sketch a proof of the symmetry ()4.19p which follows from the fact that the operators associ- 
ated to the evolutions (|4.13p . (|4.16p are adjoint. A solution f(t,x) of 

dtfit, x) = D{-p{x))dlf{t, x) + E{x)a'{p{x))d^f{t, x) 

with initial condition /(O, x), is, by definition of (j4.13p . equal to 

f{t,x) = I dyGUx,y)f{0,y). 

Now using the fact that 

f{t + dt,x) = j dyG^t+atix,y)f{0,y)= J dyGl{x,y)fidt,y), 

one can show after an integration by parts that 

dtG' = d'y[Dip{y))G'] - dy[E{y)a'{piy))G'] • 

Therefore G^{x,y) evolves according to the same equation as G\{y,x) (|4.16p and they coincide at 
time t = (Go(x,7/) = G\{y,x) = 6x=y)- Thus one concludes that they are identical. 

4.3. The case of a battery. To represent our model of a battery as a limit of the smooth 
continuum problem discussed above, we consider the field E{x) localized in a certain region of size 
Ax around the origin with E{x) = outside the region [1 — 1] U [0, ^] (where 1 is identified 
to on the ring). Then one should take the difference of the chemical potentials p{^), p{l — ^) 
at the edge of the battery to be fixed. Using the relation between the chemical potential and the 
density of systems in local equilibrium, one has 

/Ax\ ( Ax\ fP^'^^'^ 2D{p) , . , 

where the constant K' is a characteristic of the battery which remains fixed as Ax 0. 

For K' > 0, the picture is that for Ax small, the density has an abrupt increase for x in the 
region where the field is localized and a slow decay for ^<x<l — In the limit Ax — )• 0, the 
steady state profile p{x) satisfies —D{p{x))p' {x) = ^7 for < x < 1 with a jump from p~ = p{\) to 
p'^ = p{0). The current J' and the densities p~ and p"*" at the edges of the battery then satisfy 

/ -^dp = K\ / p{x)dx=-p, / D{p)dp = J. 
Jp- (^{P) Jo Jp- 

If we assume that in the limit Ax — )• 0, the Green's function Gl{x,z) converges to Gl{x,z) then 
expression (j4.20p becomes 

C(..y) = -B(.)Bfa,/..^g|||-J^~rft^^.6K-,=)4(^g||)6;fe.^) 



OO 



-J / dt 







2D{p+) 2D{p' 



Glix,0) Gliy,0) . (4.22) 
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It is then easy to check that this expression of the pair correlation functions, for a general diffusive 
system, reduces to the one of the battery model ()3.37p when one takes D{p) = 1 and a{p) = 
2pil-p). 

5. Comparison with the open system and with numerical simulations 

5.1. Comparison with the open system. For open diffusive non-equilibrium systems, the 
steady state has long range correlations which scale like 1 /N for a one-dimensional chain of length 
A^. The two-point correlation function of the SSEP in contact with reservoirs at densities pa, Pb is 
given by (SpH IDLS2] 

with 



1 ^ « < J ^ iV, iVi; ??j)opon = TrCopcn 



X<y, Copen{x,y) = -JopcnX{l-y), (5.1) 

where jTopen = Pa — Pb stands for the steady state current. 

The SSEP on a ring driven by a battery can be interpreted as the canonical non-equilibrium 
counterpart of the SSEP driven by reservoirs. By tuning the mean density p and the asymmetry 
{p,q), the boundary conditions at the battery p^,p~ (12. 9p can be chosen to be equal to Pa,Pb so 
that both systems have the same steady state density profile and the same steady state currents 
J' = jTopen- In both models local equilibrium holds and implies that at the leading order both 
steady state statistics are given locally by a product Bernoulli measure with the same density. In 
this sense, equivalence of ensembles is also satisfied in non-equilibrium. In order to understand 
more precisely the interplay between the non-equilibrium correlations, the canonical constraint and 
the driving mechanism we are going to compare the correlation functions (j3.37p and (|5.ip (as well 
as K20\\ ). 

Expression (j5.ip is related to the inverse of the Laplacian A^''' in [0, 1] with Dirichlet boundary 
conditions. Let G^^"^ = exp(— sA^^'') be the corresponding Green's function, then we know from 



Spl that 



Copcn(x,y) = -j2(AD-)-\:E,y). 



This can be rewritten as 



poo poo pi 

Copenix, y) = -2J^ / ds G^^[x, y) = -2J^ / ds dz Gf\x, z)Gf\y, z) , 
JO Jo Jo 

where J = Jopen- This is reminiscent of the expressions derived for the battery model (|3.37p 
and the weak field (j4.20p with two notable differences. The Green's functions are not associated 
to the same operator (they differ by their boundary conditions or the bulk field). In (j3.37p the 
contribution of the battery is given by the term Gl{x, 0)Gj(y, 0) which is absent in the open case. 

For a general field the current J is replaced by dz ( 2D{p(z)) ) (14.20P . 

The response to a small drive of the open and close systems is different. For a small current J , 
when the mean density p 7^ 1/2 then the long range correlation (j3.37p is of order J . Instead (j5.ip 
scales like j"^. To see this, we first note that in the battery model the mean density is linear and 
/J dxp{x) = = p. Thus 
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where the mean current is given J = — p . Expanding (|3.37p to the first order in J leads to 

C^(x,.) = -^(l + ^^(l-x-.)), (5.3) 
where we used the fact that a = 1 — '^J^ J to first order in J . We note that this coincides with the 



first order expansion of the correlations (j6.2p for the canonical measure associated to the product 
measure. At the order j'^ the non equilibrium contributions appear and (j3.37p no longer matches 
with (EH. 



5.2. The battery model at density 1/2. An explicit expression of the two-point correlations 
can be obtained by plugging in (j3.37p the expression (j3.33p of G\. For general mean density, this 
expression is rather complicated. However, several simplifications occur when the mean density is 
equal to 1/2 and in the rest of this section we will focus on this case. 

At mean density equal to 1/2, the relation (j5.2p implies that the densities at the battery satisfy 

p~ = \ — p'^ that a = = 1 s-'^cl the boundary conditions (j3.30p . (|3.32p simplify. In particular 

Gl{x,y) = Gl{x,y) and the Green's functions are associated to the Laplacian on the ring [0, 1] 



k=l 



G\{x,y) = 1 + 2^exp ( - (IkTrft) cos(27r/c(x - y)) . 

As G\{y,z) = G\{z,y), one gets 

ds dzGl{x,z)Gl{z,y)-G\{x,^)G\{y,{))= I ds GI{x , y) - G^x , 0)Gl{y , 0) 



_ -s-^ cos [2-11 k{x — y)) cos {2tt kx) + cos{2Tr ky) ^ cos(27rna;) cos(27rA;y) 

From the identity 

cos (27rkx) 1 , ^1 
VxG0,1, 2 > — V-ttt^ = — x(l-x) + — . 
^ ' ^ ^ (2kiry 2 ^ ^12 

We deduce that for a; < y 

^OO /*1 /'OO 

ds \ dzG\[x,z)Gl[z,y)-G\{x,^)G\{y,^)= / dsGl^x.y) - G\(x,^)G\{^y,^) 
Jo Jo 

1, ,11, ,1, ,1, cos(27rnx) cos(27rA:v) 

A:,n ^ 1 V / \ ; 

Finally when the mean density is equal to 1/2 then (|3.37p can be rewritten for x < y 

cos(27rnx) cos{2-Kky) 



G{x,y) 



la{p^) + ^j'-j'(l{x + y){l-{x + y)) + x)+8j' f2 



(27r)2(A;2+n2) ' 
(5 



We stress the fact that the correlation function diverges when both x, y approach the battery. 

In figure [H we compare the exact expression ()5.4p to the results of a simulation, for = 64 and 
A^ = 128 for 25 x 10^ updates per site. The agreement is very good. As our data for size A^ = 128 
lies between the data for A" = 64 and the theoretical prediction, it is reasonable to believe that the 
prediction does represent the N oo limit. 
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Clx, x] 




0.2 0.4 0.6 0.8 1.0 



Figure 1. At mean density p = 1/2, the correlation N{r]i; ?7i+i) between the sites i and i + 1 versus 
X = i/N is represented for systems of sizes 64 and 128 with 25 x 10^ updates per site. The dashed 
line represents the theoretical prediction (|5.4p of C(a::, a;). 



6. Some exactly solvable models 

6.1. Zero range process. We study now the ZRP driven by a battery defined in section [2131 If 
p = q then the dynamics is reversible wrt the product measure (^f^iV^, where the measure at each 
site is given by 

where is the normalization factor and we used the convention ^(O)! = 1. 

We are going to show that for p ^ q, the invariant measure is a product measure with a canonical 
constraint. A similar property was already derived for zero range models in contact with reservoirs 
|MF] . For any /3 such that {q — p)f3 > 0, we introduce 

= ^illZ^^,:, with ^■=p(^(i-l) + Pi^^-A±l^_ (6.1) 

The parameter /3 can be adjusted to fix the mean density. 

Proposition 6.1. The invariant measure of the zero range process on a ring of length N with 
M particles is the measure Uj3 (j6.ip conditioned to a total number of particles equal to M . The 
conditioned measure is independent of (3. 

Proof. The generator of the zero range process is given by 

Lfiri) = Ei^i,N 9im) [(/(^(^■^+^)) - firi)) + (/(r?(^'-^)) - /(r?))] 

+g{m) - fiv)) + ifir]^''''>) - fm 

+9{VN)[p{f{v^''^'^) - fm + - fiM 

where rf'^ is the modified configuration after a jump from i to j. 
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In order to check that i^^ is an invariant measure, it is enough to prove that for any function /, 
we have ^^{Lf ) = 0. 



•i - (t>i+i) 



9^ - (Pi+1 

4>i - (t>N 



aivi+i) + aivi) 1 + 



Qgim) + pgivi) i + 



5V - (pl 



- 9{m) 
pgivN 



^i+i 



+T^fi{v)f{^) g{m) [q-p + g^l— ^ ] + g{j^-^^ [p-q + p^^ 



4>N 



Since (pi varies hnearly according to (|6.ip . we obtain 



iyp{Lf) = up{f{r,) /3 



g{m) g{m) 



4>N 



+ g{m) 



■pCpN 



givi 



p4>n - q4>i 



0, 



4>N / \ (Pl 

The last equahty follows from (j6.ip which implies that qcpi —p4>N = /5- 

The number of particles being preserved by the dynamics, the invariant measure is obtained by 
conditioning z/^ to have M particles. □ 

Proposition 16.11 can be recovered from the macroscopic approach. The diffusion and conductivity 
coefficients of the ZRP satisfy D{p) = a'{p)/2. Therefore in the case of a slowly varying weak field, 
the non-equilibrium part of the correlations C^; in ()4.20p vanishes and only the canonical constraint 
remains. 

6.2. Independent variables and canonical constraints. We are going to compute the scaling 
limit of the truncated two-point correlations in product measures with canonical constraints. By 
analogy with the measures found in the previous section, we consider a chain of N independent 
variables {?7j}o < i ^ v with densities pi = p{i/N) which vary slowly from site to site. 

Let {t'pip be a family of measures with different chemical potentials which are indexed by their 
density p. The variance of the measure fp will be denoted by x{p)- We consider the measure p 
which is the product measure ^f^i^p^ conditioned to have a total number of particles equal to 
the integer part of pii / N) . The density profile of the measure p is given by p{x) and the 

two-point correlation function scales for large and x ^ y m. [0, 1] as 

1 x[p{x))x{p{y)) 



IJ'{riNx ; r]Ny) 



.2) 



These asymptotics are very different from the long range non-equilibrium correlations of the SSEP 
with a battery and the same density profile (j3.37p . 

To derive (j6.2p , one can approximate the original system in the large limit by a chain of Gauss- 
ian independent variables under the canonical constraint. Define p to be the Gaussian measure of 
N independent Gaussian variables {A"i}i s; j ^ v with variance Xi = x{p{i/^)) at site i under the 
canonical constraint "^ZiLi -^i — 0- 
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To compute the marginal of /i over the two variables {Xi,X2}, we integrate over the N — 2 



remaining Gaussian variables. Given {Xi,X2} the probability density that YliL^^i = ~ ^2 
is proportional to exp ( — y^{Xi + X2)'^) with x = Xli^s Xi- Thus the marginal of /i over the two 
variables {Xi,X2} is Gaussian with correlation matrix 

-1 

_ 1 / Xi{x + X2) -X1X2 




X + Xi+X2\ -X1X2 X2{x + Xi) 
For large N, x — X x{p{f ))dr, so that the two point correlation scales like 

f^{Xv,X2)'- ^ 



^ Jo X{p{r))dr 
Thus dlJ]) follows. 

7. Conclusion 

In this paper, we studied several NESS with a canonical constraint on the particle number and 
computed the two-point correlation functions by using a macroscopic approach. We considered two 
driving mechanisms: either a microscopic strong field localized on one bond (the battery model), or 
a weak field with an intensity which varies smoothly on the macroscopic scale. The battery model 
can be recovered as a singular limit of the smoothly varying field. Different driving mechanisms 
lead to different stationary measures with different long range correlations. In fact, the long range 
correlations of the canonical models differ also from the correlations in systems maintained out of 
equilibrium by reservoirs. 
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8. Appendix I: Green's functions 

In this appendix, we derive the exact expressions (|3.c{3|) . (|3.34p for the Green's functions G^, 
which satisfy (f330D . ([3:32]) . We recah that a = 

We start with the computation of G\ introduced in (j3.32p which evolves according to 

dtG\{x,y)=dlG\{x,y), with G\{l,y) = aGl{^,y), d,G\{l,y) = d,Gl{^,y) , (8.1) 

with GQ{x,y) = 6x=y The Laplacian with boundary conditions (18. Ih at the battery is not self- 
adjoint, but it can be decomposed on the following basis. We define 

fl{x) = {ax + l-x) cos{2kT:x), fl{x) = sin(2/c7rx) . (8.2) 

These functions satisfy the boundary conditions (18. ip . Furthermore 

dlflix) = -{2k^ffl{x) + 4A;^(1 - a)/|(x), dlfl{x) = -{2knf fl{x) . 
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Let flit, x), f'^{t, x) be the solution at time t of the heat equation with boundary conditions (|3.28|) 
and initial data fl{x),f^{x), then 



'flit, x) = exp ( - i2k7:)H) (/^(x) + 4fc7r(l - a)tfiia 



J^,it,x)=expi-i2kn)H)f, 



qix) 



(8.3) 



The Green's function introduced in (jS.SOp satisfies 

dtGfix,y) = dlGlix,y), with G?(l, y) = G?(0, y), ad,Glil,y) = d,G'^iO,y), (8.4) 

with Gg(x,y) = 6x=y For the Laplacian with boundary conditions (j8.4p . we can introduce the 
basis 

glix) = cos(2/c7rx), dki^) = ia + x — ax) sin(2A;7rx) . (8-5) 

and one has 

dlglix) = -i2kT:fglix), d^gHx) = -i2k7rfglix) + 4fc7r(l - a)glix) . 

If glit, x),glit, x) denote the solution at time t of the heat equation with boundary conditions (j8.4|) 
and initial data g\.ix),gHx), then 

{glit,x)=e^vi-i2kT^?t)glix) 

\glit,x) = exp ( - i2kT:ft) (4A:7r(l - a)t glix) + glix)) 
Using the fact that for x,y G (0,1) 

^ cos(2n7r(x -y)) = - + -S^^y . 

n=0 

one can check that 

2 ^ 4 

oo 

= ^^(«2/ + 1 — y) cos(2/c7ry) cos(2A;7rx) + (a + x — ax) sin(2/c7ry) sin(2/c7rx) . 
fc=o 

oo 



fc=0 

This implies 

-1-11 °° 

^5x=, = ^f^iy)g'oix) + J2 fliy)9li^) + fliy)9lix) . (8.7) 

fc=l 

Thus the Green's function G^ is given by 
" + l^l^„„,^ _ 1^1 



Glix,y) = -f^ix) + Y,flit,x)giiy) + flit,x)gliy) (8.8) 

fc=i 

^ oo 

2/d(^) + E^^P ( - {fli^)9liy) + fiix)gliy)+4k7Til - a)tflix)gliy)) . 



k=l 



20 



T. BODINEAU(i), B. DERRIDAt^), J. L. LEBOWITZ^^) 



Similarly by exchanging x and y in the r.h.s. of (I8.7p one can see that 
^GUx,y) = -f^y) + 5]exp ( - i2kn)H) [fl{y)gl{x) + fl{y)gl{x) + 4A;vr(l - a)t/2(y)5i(, 

A;=l 

This proves the expression ()3.33p and the identity G^{x,y) = G\{y,x) claimed in ()3.34p . 

9. Appendix II 

In this appendix, we are going to derive the relation (j3.25p from (j3.35p . 

We first start by proving some identities. The correlation between x and w is symmetric so that 
for any x and w 

dyGl{x,y)dy{a{p{y))dyGl{y,w)^ = dy G\{w,y)dy{a{p{y))dyGl{y,x) 
Using the fact that G].{x,y) = Gl{y,x) ()3.34p . one gets 

j\yG\{x,y)dy{a{p{y))dyGl{y,w)) = j\y Gl{y,w)dy{a{-p{y))dyG\{x,y)) . (9.1) 
Expanding each terms of (|9.ip leads to the two identities, 

' dy G\ {x, y)dy {a{-p{y))dyGl{y, w)) (9.2) 

dyGl{x,y)dya{p{y))dyGl{y,w) + / dyGl{x,y)a{p{y))d'^Gl{y,w), 
Jo 

dydyGl{x,y)dya{p{y))Gl{y,w) + / dy dlG\{x,y)a{p{y))Gl{y,w) . 

Jo 

We introduce the notation [F(y)]yZQ = F{1) — F{0). Integrating by parts the first term in the 
last equation of ()9.2p . one gets 











1 

dydyG\{x,y)dya{p{y))Gl{y,w) 



\yGl{x,y)dy(dya{-p{y))Gl{y,w)) + [G\{x,y)dya{-p{y))Gl{y,w)t:}^ 







\yGl{x,y)dya{p{y))dyGUy,w) - [' dy Gl{x,y)d'ya{p{y))Gl{y,w) 
Jo 

+ [Gl{x,y)dya{p{y))GUy,w)]lzl. (9.3) 

By summing the two equalities in ()9.2p and simplifying thanks to the identity (j9.3p . we finally 
obtain 

2 / dy Gl{x,y)dy(a{p{y))dyGl{y,w] 







\ydlG\{x,y)a{p{y))Gl{y,w) + [\y Gl{x,y)a{p{y))d^yGl{y,w) 
Jo 

^ dyG]{x,y)dla{-p{y))Gl{y,w) + [Gi(x, y)a,a(p(y))G^(y, t.)]g . (9.4) 
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Note that dya{p{y)) = — 4j7^^, where J7 = p"*" — p~ is the macroscopic mean current (for large 
i7 = Nq with q introduced in (12 .Gp ). The discontinuity of the density at the battery imphes 
dya{p^) — dya{p'^) = — 4(p+ — p'')J = —Aj'^. Thus integrating (|9.4|) wrt time leads to 

ds j\yG]{x,y)dy{a{-p{y))dyG'l{y,w)) (9.5) 



JO 



\yG\{x,y)'^^P^Gl{y,w) - \a{-p[x))5^=^ 

+2J^ r ds f\yG]{x,y)Gl{y,w)-2J^ ds G\{x,^)G''Mw) . 
Jo Jo Jo 

Letting T go to infinity, one gets from the expression of G^ (j3.33p and G^ (13.34p 

lim [\yGUx,y)^^^GUy,w) = -^-^a{x)a{w) [\ya{p{y)), 
T^oo Jq 2 (a + ly Jo 

with the notation a = ^|^Ty a{x) = {ax + 1 — x). 

As the steady state is linear, a{p{y)) = o-{p'^) + y{a{p~) — cr{p'^) + 2J''^) — 2j'^y'^, and we get 

hm / dyG\{x,y)^i^G'T{y.w) = a[x)a{w) ' + ^ • 

Using (19. 5j) and the symmetry G\{w,z) = Gl{z,w) (I3.34p . we have finally shown that 

lim I ds I dyGl{x,y)dy(a{p{y))dyGl{y,w)) = -l;a{p{x))6^=uj - C{x,w) , 
^°°Jo Jo ' ^ J I 



T^oo 

with 



C{x,w) = -^-l^(^^^Pll±^i^ + !^yax + l-x)iaw + l-w) (9.6) 

roc rl 

-2jM ds dzGl{x,z)Gl{w,z) -Gl{x,0)Gl{w,0) . 



The previous expression combined to (I3.35P leads to ()3.36p . 

10. Appendix III 

10.1. Uniqueness of the steady state. We will show the uniqueness of the smooth solutions of 
(031) 

d,(^D{p{x))p'{x)-E{x)a{p{x))) =0, (10.1) 

for regular coefficients and D positive. 

Suppose that pi{x), p2{x) are two solutions of (jlO.ip with the same mean density. Then there is 
i7i > J2 such that 

-D{pi{x))p[{x) + E{x)a{pi{x)) = Ji 
-D{p2{x))p'^{x) + E{x)a[p2{x)) = J2 
If the profiles coincide at xq then substracting both equations, one has 

L»(pi(xo)) [A{xo) - p'^ixo)] = J2-Jl<Q- 

This implies that when the two solutions cross then pi > p2 before the crossing and pi < p2 
after the crossing, so that the solutions cannot cross more than once. But the profiles are smooth. 
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periodic and have the same mean density, thus they have to cross an even number of times. This 
is a contradiction and therefore both profiles have to be equal {^J\ = J2). 



10.2. Convergence of the linearized evolution. We turn now to the convergence of the lin- 
earized evolution (j4.16p 

X) = dl [d{p{x)) fit, x)) - {E{x)a'{-p{x))f{t, x)) , (10.2) 

with smooth initial data /(O, x) and mean dx f{0,x) = 1. It is equivalent to consider the 
evolution 

dtfit, x) = {d{p{x)) dj{t, x)) - (a{x)f{t, x)) , (10.3) 

with a{x) = —dxD{p{x)) + E{x)a' {p{x)). At any time t > 0, f{t,x) can then be interpreted as 
the probability density of a particle evolving on the ring [0, 1] with a non homogeneous diffusion 
coefficient D{p{x)) and a drift a{x). As the particle evolves on a compact set, it will reach a 
stationary state when t goes to infinity. The limiting density for the particle position will be 
denoted by f{x) and it is the unique solution with mean density dx f{x) = 1 of 

dx (^D{p{x)) dj{x)) - (a(x)/») = . (10.4) 
A way to understand the relaxation to the stationary state is to check that the relative entropy 

S{t) = - [\x fit, x) log (^^^) , (10.5) 







is a Lyapunov function for the evolution (jlO.Sp . This is a general fact for Markov processes. Taking 
the time derivative one has 

dtSit) = dx [Dipix)) d^fit, x) - aix)fit, x)) log • (10-6) 



We now note that 



\xaix)fit,x)dxlog(^^;^) = C dxaix)fix)d, ^ ^^^'"^ 







fix) J Jo ^ - ^ ^ - V fix] 
dx Dipix)) dJix)d, 



fix] 

where we used (jl0.4p in the last equality. Combined with (|10.6p . this leads to 

dtSit) = ^'dxi)(Mx)) (-^a./(t,x)-a./>))a,.(:^) 

dxDipix))fix) (^dx^l^^ . (10.7) 

Thus Sit) is a Lyapunov function. A quantitative estimate of the approach to equilibrium could 
then be obtained by using a log-sobolev inequality. 
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